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Abstract 

(— ! \ 
'Tj \ We study the growth of a crystal in presence of a correlated disorder on the 

O ' substrate. Using functional renormalization group, we show, for a long range 

disorder correlation, an initial decay of the KPZ type nonlinearity, though 

- over a large length scale the behavior can be governed by the nonlinearity. 

-)— > 

t^ ■ I. INTRODUCTION 

fi 

The equilibrium shape of a crystal surface undergoes a roughening transition from a rough 
Pi ■ high temperature phase to a smooth low temperature phase as the temperature is decreased 

R ' llh0- Above the roughening transition T^, the height fluctuations grow logarithmically 

with the dimension of the system L and below the transition the height is smooth and is 
independent of L. The non-equilibrium counterpart, the growth mechanism of such crystals, 
^ ■ provides more insights about the roughening transition 0. It is found that for T > T^. the 

growth is non-activated in nature whereas for T < Tr, the growth is essentially by nucleation 



Ti 






O . of droplets and the growth velocity is exponentially slow in the inverse of the force. For an 



O 



infinitesimal force F, the mobility which is the ratio of the growth velocity and F vanishes 



t^^ I with a jump from a finite value at the transition. With a finite force the transition is blurred. 

^ ■ The flat phase is destroyed over large length scale and eventually becomes rough. In this 

c^ ■ case a continuous decrease of the mobility is observed as the temperature is lowered [^ . 

An interesting development in this direction is the inclusion of the disorder of the under- 
__i , lying substrate [Q . The morphology of the growing surface shows a new transition which is 

C3 I the super roughening transition at temperature T^^. For T > Tgr, the height fluctuation is 

^ ■ the same as the high temperature phase of the usual roughening transition but for T < T^^ 

the surface is rougher than the thermal phase. There is a controversy regarding the rough- 
. ness of this super rough phase |]^-0 though a recent numerical treatment |Q favors the size 

r^ ! dependence of the roughness as (LogLy. A dynamical renormalization group treatment |^ 

a I shows that for T > T^^ the effect of the disorder essentially vanishes over large length scales 

and the scaling properties of the surface are the same as that of a surface growing on a 
pure substrate, in the rough phase. The linear response mobility in this super roughening 
transition vanishes continuously, unlike a jump in the case of the roughening transition, as 
Tsr is approached from the high temperature phase. For T < Tgr there is a temperature 
dependent dynamical exponent and a nonlinear response of the system to an external force. 
The connection between these roughening and super roughening transitions has been eluci- 
dated furthermore by a renormalization group treatment of the growth on a substrate with 



a correlated disorder [^]. The roughening transition turns into a super roughening transition 
as the correlation of the disorder decays sufficiently fast. By tuning the power law of the 
long range correlation of the disorder it is possible to go from the pure to the short range 
disorder limit. 

Such growth problems for crystals need to take into account the periodicity perpendicular 
to the crystal surface and is usually incorporated in the growth equation by a periodic 
pinning potential. The low temperature, flat phase in the roughening transition is due to 
the relevance of this pinning potential. 

The nonequilibrium growth problem is complicated. One has to take into account a 



relevant nonlinear term that appears due to the lateral growth of an oblique surface |]T0 
One of the major consequences of this Kardar-Parisi-Zhang (KPZ) nonlinearity is a power 
law (not logarithmic) growth of height fluctuations. In the context of the nonequilibrium 
situation of a crystal surface, it has been argued that with a finite force such a nonlinear term 
is generated by the interplay of the force and this pinning potential ^^. The asymptotic 



behavior of the surface might then be determined by the nonlinearity which can destroy the 
roughening transition. A rough phase can appear even at low temperature. If one starts 
with a nonlinear term, the lateral growth and the pinning potential combinedly lead to a 
phase factor in the pinning potential [|T^. This phase factor is also renormalized as one 
looks at the system over larger length scales. The renormalization of the phase factor makes 
the identification of the phases a nontrivial problem. 

The disorder substrate case is also not well understood. Previous analysis in this direc- 
tion with a short range correlation of the disorder of the lattice showed that in the presence 
of the nonlinearity a small driving force is relevant and the asymptotic properties are es- 
sentially governed by the nonlinear term since the finite velocity leads to a smearing of the 
pinning potential |]I2|. A numerical investigation on the other hand showed a generation of 



a quenched random mobility [|T3| and hence a new universality class. The latter has been 
attributed essentially to the effect of a finite lattice cutoff. 

In this paper, we address the question of a large scale description of the dynamics of 
a growing surface in the presence of a disorder on the underlying lattice and a KPZ type 
nonlinearity discussed before. We consider a very general form of the disorder correlation. 
In view of the confficting scenario mentioned above our analysis is based on a renormal- 
ization group technique that involves a finite lattice cutoff. A standard way of describing 
such a growth phenomena is to start from the equation of motion for the height 0(r, t) at 
time t and at coordinate r in the two dimensional plane. As mentioned before the lattice 
structure orthogonal to the crystal surface is respected through the inclusion of a periodic 
potential. The equation of motion we need to study is therefore the usual Sine-Gordon 
equation |^,|9|JTTHr3| subjected to a constant driving force F 



m-^^^^ = KV^<P{r, t)-V sin[0(r, t) + d{r)] + A/2(V</.)2 + F + R{r, t) (1) 

where K denotes the stiffness of the surface, m is the microscopic mobility, V is the strength 
of the pinning potential. Here A represents the strength of the nonlinearity that allows a 
lateral growth of the surface and R represents the thermal noise at temperature T with a 
short range correlation given by 

{R{r, t)R{r', t')) = (2T/m) 6{r - r')6{t - t') (2) 



In the presence of a random substrate we expect that the minima of the periodic potential 
will be randomly shifted and this is incorporated in this equation through a quenched random 
variable d{r). We consider here a general form of the correlation associated with this random 
phase shift d{r), as gQe^W'^)-<^(''')) = 7(r — r'), where Qq = }^iv?V'^ and the overbar denotes 
average over the disorder. At this point we do not specify the functional form of 7(r — r') 
but later we concentrate on 7(r) ~ r~^". We obtain back the short range disorder case for 
large a and the perfect crystal for a = 0. 

Various limiting forms of this equations which have been studied are as follows, (i) 
The most well studied limit is ^ = [Q. (ii) The pure equilibrium growth that shows 
a roughening transition is obtained with A = F = d = |0J^. The width of the growing 
surface is conventionally described by the following scaling w{L,t) = L^f{t/L^), where 
X is the roughening exponent and z is the dynamical exponent. Above the roughening 
transition, the width scales as w{L,t) ~ ln[Lf{t/L^)], which implies x = 0. In the low 
temperature fiat phase w{L,t) is independent of L. (iii) The nonequilibrium situation of 
this pure problem is described by equation ^ with d = 0, with or without F ||lT|. (iv) 
The short range correlated disorder on the substrate follows from 7(r) = go6{r) [Q. This 
shows the continuous super roughening transition, (v) A long range correlated disorder 
with 7(r) ~ 7--2a a^^^ A = shows the interpolation between the roughening and super 
roughening transition 0. Another interesting limit, V^ = A = and a quenched noise -R(0, r) 
was used to study the dynamics of driven interface in a disordered medium [^ . 

In a renormalization group (RG) approach the system is looked at longer length scales 
by integrating out the effects of small scale fluctuations. When the system is rescaled to the 
original scale, the effect of the small scale fluctuations goes into the renormalization of the 
various parameters of the problem. The effective coupling constant observed at a certain 
length scale is then given by the RG recursion relations. In this approach, in all the above 
cases except (ii) with F y^ 0, X has been found not to be renormalized. The important 
new feature that arises in our case is the renormalization of the nonlinearity due to the 
nonlocal property of the correlator. We obtain a surprising result that this nonlinearity 
decays initially with the length scale for F = 0. However the possible generation of a force 
and its relevance can cause a truncation of this decay of the nonlinearity and the growth 
asymptotically becomes KPZ like. In order to treat the long range correlation in general, 
we use the functional renormalization group (FRG) approach where the renormalization of 
the correlation function 7(r) and its effect on the other parameters are studied. 

The paper is organized as follows. Section II is devoted to the description of the effective 
generating functional. In section III we derive the functional renormalization group flow 
equation with necessary diagrams. At this point we might add that our approach is simpler 



than the dimensional regularization approach of Ref ||T^. Details of the calculations are 
presented in the appendices. Section IV is devoted to the discussion of the asymptotic 
behavior of the system and connection with previous predictions. In section V, we summarize 
our results. 



II. GENERATING FUNCTIONAL 

The growth on a crystal substrate with disorder is described by the phase disordered 
Sine-Gordon model with a KPZ nonlinearity as described in equation |I[ We use the Martin 



Siggia Rose formalism [|T6| that requires a response field 0(r, t). Averaging over the disorder 



J7| yields the generating functional Z = fVcjyVcf) exp[^]. Here A = A^ + Aq is the 
effective action with free and the disorder part of the action given respectively by 

4°) = J dtdr{^Mo{r,tf - 0o(r,t)KVo(r,t) - «:oVVo(r,t) - ^(V0o(r, t))'] 
+Jo(r,t)0o(r,t)} (3a) 

4"^ = J dtdt'drdA-foir - rO<^o(r, t)Mr', t')cos[Mr, t) - M^', t')] (3b) 

where (po = (j), '&o = 2mT, fiQ = 1, kq = mK, Aq = Xm and Jo = niF are the bare quantities. 
The Gaussian part of the action (Eq. ^ with Aq = ) gives rise to the following response 
and the correlation functions @] in the momentum and frequency representation, 

(0(q, a;)0(q', uj')) = R{q, uj)5{q + c^)6{uj + uj'), 

(0(q, ^)0(q', cu')) = C{q, u)5{q + q')5(^ + cu'), (4) 



where, 



RiQ,^) = — — : 2^ 

c(?,^) = ^- — Y^ri (^) 



(suppressing subscript 0). In the momentum and time representation 

(0(q,t)0(-q,t'))=^(t-i')/ie-'^'^''(*-*'), 



(0(q,t)0(-q,t')) = ^e-'^'='^'l*-*'l, (6) 



where d{t) = 1 if t > and zero otherwise. In real space 



i?(r,t) = ^il^e-'^'/(^^'^*). (7) 
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The correlation function has a divergence due to long wave length and for our purpose we 
define a difference correlation with a suitable ultraviolet regularization introduced by a cutoff 
A as following. 



1 . , , 2\ H'd [ dk 



-fj.Kk^\t-t'\] 



r(| r-r' \,t-t') = -([0(r,t)-0(r',t')f) = ^ / -[1 - Jo(A: | r - r' |)e 

(8) 

Our analysis involves this difference correlation. The correlation and the response functions 
are also connected by the Fluctuation Dissipation theorem (FDT) given as 

e{t>0)dtC{k,t) = -^R{k,t) (9) 



An important aspect of our system is the breaking of the FDT. This violation of FDT has 
previously been argued in the dynamics of random phase sine Gordon model |18|. The 



diagrams corresponding to the correlation and response functions are presented in Fig 1 
along with the other two vertices (j)(f)'cos{(f) — 0') and 0(V0)^, where 0' and 0' denote the 
fields at (r',t'). These are the basic diagrams needed for the subsequent perturbative FRG. 



(a) 



AAA 
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R 



(b) 



(c) 

Figl 



(d) 



Fig 1. Diagrammatic representations of (a) correlation function, (b) response function, (c) 

0(r, t)0(r', t')cos[0(r, t) — 0(r', t')] where dots represent many lines, (d) 0(r, t)(V0)^ 



III. RENORMALIZATION 

In the following we discuss the renormalization scheme to obtain the physics at large dis- 



tances, which, in the momentum space, corresponds to small k ||T9[. We use the description 
of the fields as sum of the fast and slow modes defined as, 



(ps = 0(A;) for < A; < A, 
0(A;) for A < k < A + 6A. 



(10) 



This separation of fast and slow modes is used to average out the fast fourier modes or the 
short wavelength details of the problem. This short wavelength properties are incorporated 
by the appropriate renormalization of the parameters in the effective action which describes 
the large length scale properties. By separation of the fast and slow Fourier modes, in 
general, one arrives at the action. 



Z= I Y[ V(j)V(j)e 

0<k<A 



■A.(i{4's,<j>s) 



A<k<A+5A 



Ao{<f'f,'i>f)pAi{(t>s,<t>s,(j>f,<t>f) 



(11) 



where ^o is the free action and Ai is the interaction part that contains both fast and slow 
modes. Here 



^A'{(ps,<Ps) _ ^Ao{(l>s,<i>s) /^Ali(t>s,<l>s,<l}f,<ff)\ 



0>, 



(12) 



where < >o> denotes average with respect to fast modes of the free action. The next step is 
a cumulant expansion to clearly identify the contribution from averaging of the fast modes, 
as follows 



^A'{<t>s) ^ ^Ao+M 



(13) 



where 



5A={Ai) + {{A])-{Aif)/2 + ... 



The system in its original length scale is retrieved back again by rescaling the 
fields. Rescaling implies that under the transformation x — > 6x, t ^ 6^t, (f) — »• 
b^cf) and (j) —>■ h^cf). Around the Gaussian fixed point we have z = 2 and 
in a 2-d system the rescaling leads to — i> 6~^0 with remaining invariant. 




I 



(a) 




(b) 



(c) 




Fig 2 



Fig. 2 (a) New vertex 0(r, t)0(r', t')(V0(r, t))2cos[0(r, t) - 0(r',t')]. Hatched line 

represents V0(r,t), (b) and (c) Origin of this new vertex. Figures on the left side of the 

arrows show the two loop diagrams. The new vertex is generated by snipping off one 

internal line of the two loop diagram. Diagrams found by this procedure are shown on the 

right side of the arrows. 



The terms that contribute to the renormalization of various vertices have diagrammatic 
representations that are presented as we proceed. In the following we consider one loop 
diagrams. In fact in this approach since each internal line in the diagram is within the shell 
A + (5A it is sufficient to consider the diagrams with only one loop |20|. The effect of the 



higher loops in the field theoretic approach |T^ are taken care of in our approach by the 
new vertex generated under renormalization. This feeds back to the renormalization of the 
original vertex we start with. 
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The nonlinear term causes a renormalization of the tension k and the temperature ■^. 

We shall not elaborate on this part since they are very well documented in previous work 

I0| . The mobility fi, as well as •& and k are renormalized by the pinning potential. For local 

correlator this can be found in references [^3 • The nonlocality of the correlator requires a 



more general treatment presented in the appendix. There is a possible new renormalization 
of the pinning potential due to the nonlinearity. A straightforward diagrammatic expansion 
would produce two loop diagrams [0 as shown on the left sides of Fig 2b and Fig2c.To 
include the effect of these in our RG scheme, we have to involve new vertices whose one loop 
contribution would be equivalent to the effect of original diagrams. This new vertex shown 
in Fig 2a is generated by the combination of A and 7(r) and will be denoted in the following 
as 7i(r). The explicit form of this vertex is 



dr dr' dt dt' 71 (r - r')0(r, t)0(r', t')cos(0(r, t) - 0(r', t'))(V0(r, tjf 



(14) 



It is possible to choose a symmetrized form for this vertex but the results will remain the 
same. 

The nonlinearity also gets renormalized unlike the case with a local correlator of Ref 
1^ . The diagrams contributing to the renormalization of the nonlinearity is shown in Fig.3. 




(a) 




Fig. 3 Diagrams contributing to the renormalization of A. 



The details of the derivation of various terms present in the recursion relations are described 
in the appendix. In appendix A, we show the necessary details to find out the renormalization 
of the tension. In appendix B and C, the renormalizations of the disorder correlator and the 
nonlinearity are discussed respectively. The final flow equations of the parameters in terms 
of a dimensionless quantity rii 



■^ are given as 



Ok 



inn 



drA{r) 



(15a) 



~dl 



4:ni 



dt drB{r,t) 



(15b) 



dd ^ f . r r^, ^ ^V^ 



^n.JdrdtBir,t) + ^J-- (15c) 



dl J ' ' ' 16K3(27r) 

= (4 + r ■ V)7(r) - 2ni{-f{r) - 2A'^-fi{r)) (15d) 

(2 + r • V)7i(r) - ^n^(7(^) - 47i(r)A') (15e) 

drA{r) (15f) 



^^ = (2 + r.V)7i(r)-^^^, 
aA A 



dl IQtxk^ 

dJ 

— = 2J + ni\ (15g) 

where A{r) = 7(r)r^Jo(Ar)e~^*^'''°) and B{r,t) = 7(r) Jo(Ar)e~'^''^ tg-y(r,i)_ r^^^e flow equa- 
tion for 7 in tliis form apparently differs from that of ref. [jl2| for a short range correlation. 
They obtain a A dependent contribution which in our formulation appears through a new 
vertex 71 generated under renormalization. Since 71 = is the initial condition, one can 
write for simplicity the flow of 7 as 

^^ = (4 + r ■ V)7(r) - 2n,7(r) - CA27(r), (16) 



where C is a constant as it is in ref. ||T2| 



IV. DISCUSSION 

For A = 0, the growth equation, Eq. |l], corresponds to the equilibrium problem 0. In 
this case, we recover the equations of Ref. . The recursion relation for the correlator of the 
form 7(r) ~ Cor~^" for large r, implies that a remains unchanged under renormalization, 
as has been observed in ref. 0. It is only the amplitude Cq that is renormalized. It is also 
apparent from Eq |15d| , that for A = the relevance or irrelevance of 7(r) is determined by 
the value of rii, that plays the role of temperature. For a perfect crystal with 7(r) = qq, 
rii = 2 determines the roughening transition whereas for short range disorder, 7(r) = gQ6{r), 
this boundary is at n; = 1 0,^. 

For 7(r) = goS{r) the flow equation shows that A is invariant under renormalization. This 
is in agreement with ref |12|. For any other short range form of 7(r) (e.g. exponentially 
decaying 7(r)), it is possible to see from Eq 0, that the width of the correlation defined as 
J drr'^'-f{r)/ J 'y[r)dr vanishes as one approaches I -^ 00. Interestingly, for a long range form 
of 7(r) as mentioned above the integral converges for a > 1. This a = 1 is the borderline 
for the roughening and super roughening transition |Q. We therefore conclude that any 
short range form of 7(r) leads to the same conclusion as ref. [|T2|, namely that the large 
scale behavior of the growing surface is KPZ like and the super roughening transition is lost. 
With a power law correlation of the disorder, this short range limit is achieved when a > 1. 
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Let us now consider < a < 1. Note that the integral over r in the flow equation 
for A is the same as that in the flow equation for k. Therefore with 7(r) ~ CQr~'^°', this 
integral is positive as long as n > 5/4 — a as in ref 0. As a consequence of this the 
nonlinearity decreases with the length scale. This is quite an astonishing feature arising 
due to the nonlocality of the correlator. Recall that the A term occurs to take into account 
the lateral growth of the surface. Our result suggests that at least at the initial stage, 
the lateral growth of the evolving surface is inhibited by the long distance correlation of 
the quenched substrate disorder. The effect of this decrease of A implies that the super 
roughening transition is unaffected up to the length scale below which our treatment with 
a F — > is valid. If we stop at this leading order of the flow equations, then from 0, we 
can argue that for a < 1 one would observe roughening transition. There is, however, a 
possibility of generation of a driving force due to the nonlinearity A as given in the recursion 
relation |15g| . The growth of this force with the length scale requires a non perturbative 
treatment of the force []TT| (explained below) and the decay of the nonlinearity is eventually 
prevented. The relevant length scale at which this crossover takes place, can be obtained by 
solving the coupled differential equation with and without the force term. Also note that, 
to this order, the flow of A and k, under renormalization, has an invariance Xn'i = constant, 
the physical origin of which is not clear to us. 



Let us consider a finite force case. According to reference |11| a finite force for a perfect 
crystal (a = 0) implies a rough surface with the height fiuctuations having a power law 
scaling with the length. To treat a finite force we redefine the height as 

0(r,t)^0(r,t) + Ft 

where F is the external force. This implies a replacement of the disorder part of the action 



as 



Ai^^ = J dtdt'drdA^oir - r')0o(r, t)0o(r', t') cos[0o(r, t) - 0o(r', t') + (t - t')F] (17) 

The next step is the splitting of cosine into a product of sine and cosine terms. Most 
crucial is the renormalization of A which follows from the possible contraction of the term 
sin[0o(r, t) — (/>o(r', t')] sin[(t — t')F]. This leads to the fiow equation for A analogous to that 
in the reference |TT[]. We do not go into further details of the calculation since the procedure 
for extracting this term is similar to the details given in the appendices. 

The connection with reference |]T3[ can also be understood in our approach. There the 
equation of motion contains a random force. Generating functional obtained after averaging 
over this random force is in fact the first term of Eqn. ^ in the expansion of the periodic 



cosine function. The random lateral drift velocity, as suggested in ref. [13|, appears through 
a new vertex 71 of Eq. |l^ again as a first term in the expansion of the periodic function. 
The random part of the KPZ nonlinearity is not present in our work but it will obviously 
be generated in higher order of A contributions. In other words, the terms that have been 
discussed in Ref. [|13| are naturally generated or present as the leading order terms in our 
functional renormalization scheme with a finite cut off, though they do not occur in the 
dimensional regularization scheme of Ref. [Q . We have treated the periodic function com- 



pletely. It is clear from section III that our renormalization is crucially dependent on the 
periodic function, and in fact, truncation to first order of the cosine will never lead to any 

9 



renorinalization of A. Similarly, the finite force scheme is also dependent on the periodic 
function. We are therefore not sure whether the leading order terms in the expansion of the 
cosine function of Eqs. pE and IH can capture the whole effect, especially the question of a 



new universality class. To address this question one has to go to higher orders in RG. This 
remains to be done. 



V. SUMMARY 

In this paper we investigated the super roughening transition in presence of a nonlinearity 
supporting a lateral growth of the surface. The disorder is considered to have a long range 
correlation as discussed after Eq|^. A functional renormalization group analysis using a finite 
momentum cutoff leads to flow equations for the parameters of the growing surface described 
by Eq. |I| and ^, and a flow equation for the disorder correlation. These equations ( [L5a[ - |15gj) 



predict the macroscopic properties of the surface. Our equations reduce to the known form 



of reference [0] for short range correlation of the disorder. The lateral growth governed by 
the nonlinear term with coefficient A is initially suppressed due to the long range disorder 
correlation. This is apparent from equation |151| . Therefore at least at the initial stage one 
might expect to see that the super roughening transition remains unaffected. Over a large 
length scale a finite force can ultimately lead to a rough surface and destroy the transition. 
This requires a treatment of a finite force as discussed in section IV. In our analysis, we have 
taken a finite cutoff and have shown that the terms that have been predicted from numerical 
lattice simulations are really generated though our analysis envisages a more general term. 
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APPENDIX A: RENORMALIZATION OF THE TENSION 

We discuss only the contribution of the pinning potential 7(r) to the renormalization of 
the tension n. The renormalization of the other quantities like temperature, and mobility 
can be found out in an analogous way with appropriate diagrams. One representative Feyn- 
mann diagram that contributes to the renormalization of the tension is the shown in Fig 4. 




Fig. 4 Diagram contributing to the renormalization of the tension k. 
The external lines are associated with momentum K and all the internal lines contain fast 
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momenta. However the sum of all of the internal momenta should add up to yield K. Let 
us choose one internal C line with momentum A as shown in Fig 4. Its contribution in real 
space is 

jdp A(5(p- A)exp(zp-r)C(p,t) = ^Jo(Ar)^exp(-/x«:A2 | t \) ^ G*"^(r,t). (Al) 

Therefore the contribution of Fig4 is 

- j dt f dr f dkexp[-i{K^k)-r]G^'^{r,t)R{r,t)exp[-Y{r,t)]-f{k) (A2) 

where exp[— y(r, t)] takes care of the sum over all possible diagrams that arise from all 
possible contractions of the internal lines as it is in the case of local correlator [|,|l3]. The 
other possibility arises when one associates the R line with the momentum A. In this case 
we have the line with the fast momentum as 

/ dp AS{p - A) exp(zp ■ r)Ro{p, t) = — Jo(Ar)/ie-'^'^^'* = C^^^ir, t), (A3) 

and the full expression is 

- fdt /"dr /"dkexpH(K-k)-r]G"^^(r,t)exp[-r(r,t)]7(A;) (A4) 



The flow equation for the tension is obtained by addition of Eqs. |A4| and |A2| with an extra 
symmetry factor, completing the integral over k and extracting the term proportional to 
K"^. Simplified form given in Eq. p.5a| is achieved by the use of the FDT and a subsequent 
integration by parts over t. 



APPENDIX B: RENORMALIZATION OF 7 

The first term in the flow equation of 7 is, in fact, the result of the contribution that 
comes from expanding the cosine in the pinning potential. The simplest diagram in the 
one loop level is shown in Fig 5. 



/\AAA^ 



A/W 



Fig. 5 Diagram contributing to the renormalization of 7 



The corresponding expression is 

7(r) f dk 



Jj^,A6{k-A)C{k,0). (Bl) 



This leads to the first term of the RG contribution in Eq. |15d[ The second term with 71 



in Eq. |15d] will be clear from the following discussion. The contribution that comes from 



the combined effect of the nonlinearity is at the two loop level and the possible diagrams 
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are shown explicitly in Ref [Q. The effect of these diagrams are incorporated through a 
higher order vertex. Usual procedure that is adopted is to snip off the internal line [^ and 



generate a new vertex. This vertex will renormalize the vertex 7 when they are rejoined 
again. We pick up one such two loop diagram which contributes in the field theoretic 



approach. Snipping off one internal line leads to a vertex of the form given in Eq. |T^. 

We call this new vertex 71 (r) which differs from the 7(r) vertex by two V0 type legs. 
Once they are rejoined this vertex will contribute to renormalization of 7(r). The second 
term with 71 in Eq. [15d| follows from this prescription. We now discuss the renormalization 
of the vertex 71 . There are two such possible diagrams which can originate from the original 
two loop diagrams and lead to the renormalization of 71. From the other two loop diagrams 
it is easy to observe that no such vertex that respects the symmetry of the system can be 
formed. The contribution of the 71 vertex in the renormalization of 7 is identical to the 
contribution of 7 vertex itself except for a multiplicative factor A^ which is obvious from a 
dimensional analysis. As it is evident from Fig 2, this new vertex is generated from the term 
of 0(A^7). It is easy to evaluate the expressions of the diagrams in Fig2. In the momentum 
and frequency representation, 

\^ roo j-t2 r 

Fig2b = -—-fir) J dt2 J dti J rfA;2^(/^2mi, ^2 - ti)C{k2mi, t2)R{k2pi, h) 
(ki ■k2pi)(ki -ksmi) 

'-kl (B2) 



/^^ A^ 7(r) ^^2 



where k2pi = k2 + ki/2, and k2mi = k2 — ki/2 in the hydrodynamic limit i.e fci -^ 0. 
Similarly 

F.,2c=^^^4 (B3) 

A similar contribution comes from the vertex involving 71. The combined expression leads 
to the recursion relation for 71 in Eq. p.5e| . The new vertex 71 is completely the effect of 
renormalization. 



APPENDIX C: RENORMALIZATION OF A 

We notice that the possible one loop diagrams which contribute to the renormalization 
of the nonlinearity A in the case of nonlocal correlator are those shown in Fig 3. There 
are other possible diagrams that could lead to a contribution in the renormalization of the 
nonlinearity. But we note that they cancel altogether. For example the diagrams shown in 
Fig 6. The expression corresponding to the diagrams in Figure 3a is 

FigSa = -^ / —^duiR{ki,Ui)R(ki, -Ui) / dr / dt^{r) ^g-^^'-'^^e^'^^-'-e-'"^*. (CI) 
4 J (27r)^ J Jo 2 

After the integration over the shell momentum ki we obtain 
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Furthermore the integration over the frequency leads to the final contribution of Fig 3a as 



dt dv r^ 7(r) Jo(Ar)e^^("'*)e-^''^ *. 



?)1'KK Jo 

Similarly the expression corresponding to Fig 3b is 



(C2) 



Fig3b = 



(iki 



Fig 3b = - 
An use of the relation 



duJiR{ki,uJi)C{ki,uJi) dr dtr^R{r,t)-f{r)kye'''''e~^^''''\ (C3) 



(C4) 



16 J (27r)2 J JO 

After performing the integration over the momentum shell and the frequency, we get, 



A^/i 



1287r/t2 



dr 



dtr^-f{r)R{r,t)Jo{Ar)e 



-fiKA'^t -Y{r,t) 



d{t > 0)dte 



-Y{r,t) 



-^i?(r,t)e-^(-*). 



(C5) 



and subsequent integration over t leads to the following expression 

A 



+ 



Fig 3b 

XfiA^ 



QAtik? 



drr^Jo{Ar)-f{r)e-^^'''^'> 



dr I rftrVo(Ar)7(r)e-^("'*)e-'^''^'*. 



, , „_ , _._. _,._.,,,.,_ _ (C6) 

GAttk J Jo 

Combining the contributions of Fig 3a and Fig 3b with appropriate symmetry factor we 
arrive at the recursion relation for A. 

The diagrams corresponding to Fig 6 can be evaluated in a similar manner. 



Vvv-s—Y. . -Kap 




(al) 



(a2) 





(b2) 



Fig. 6 Diagrams which could contribute to the renormalization of A but vanish due to their 

mutual cancellation. 
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FigGal = J — -^ j dujR{p+, uj)C{p^, uj)p ■ p_ 

/•oo 

dr dt-f{r)R{r,t)Jo{p-r)e-^^''^'^e'^\ (C7) 

Jo 



—-- j dujR{p+, uj)C{p-, cj)p ■ p_ 

/•oo 

dr rft7(r)i?(r,t)e-*^*e~^('^'*)Jo(p+r) (CJ 



where p+ = p/2 + k, and p_ = p/2 — k. In the above expressions p is the slow momentum 
and the contribution to the A vertex can arise from the following terms after adding Fig 6al 
and 6b 1 

JdvJ^^dt 7(r)Jo(Ar) R{r,t)e'''^^^'^ J ^2p^ p ■ (p/2 - k) 



du e 



iujt 



{pi + UJ^){pt + UJ^) 

+ JdrJ^^dtr 7(r) Jo(Ar) R{r, t)e-^(^'*) J ^p ■ (p/2 - k) 

{p% + UJ^){pt +UJ^) 

We need not go into the final expression corresponding to the combined contribution of Fig 
6al and Fig 6bl. This is because, interestingly, the other two diagrams in Fig 6a2 and 6b2 
lead to the same contribution but with opposite sign, thereby canceling the contribution of 
Figs 6al and 6bl. This cancellation leads to a considerable simplification of the recursion 
relations. 
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